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ABSTRACT 



u 

Context. High energy electrons can produce Hard X-Ray (HXR) emission in galaxy clusters via Inverse Compton Scattering 
Q ' (ICS) of CMB photons. 

Aims. We discuss here the consequences of the presence of such high energy particles for the multi-frequency emissivity of the 
C/3 ^ same clusters and for the structure of their atmospheres. 

. Methods. We derive predictions for the ICS HXR emission in the specific case of the Ophiuchus cluster under three main scenarios 
for producing high-E electrons: primary cosmic ray model, secondary cosmic rays model and neutralino DM annihilation scenario. 
We further discuss the predictions of the Warming Ray model for the cluster atmosphere. Under the assumption to fit the HXR 
^ emission observed in Ophiuchus, we explore the consequences that these electron populations induce on the cluster atmosphere. 
Results. We find that: i) primary electrons can be marginally consistent with the available data provided that the electron 
spectrum is cutoff at £ i$ 30 (90) MeV for electron spectral index values of 3.5 (4.4); ii) secondary electron models from pp 
collisions are inconsistent with the viable gamma-ray limits, cosmic ray protons produce too much heating of the intra cluster 
(IC) gas and their pressure at the cluster center largely exceeds the thermal one; iii) secondary electron models from DM 
. annihilation are inconsistent with the gamma-ray and radio limits, and produce too much heating of the IC gas at the cluster 
' center, unless the neutralino annihilation cross section is much lower than the proposed value. In that case, however, such models 
^\ _ no longer reproduce the HXR excess in Ophiuchus. 

r Conclusions. We conclude that ICS by secondary electrons from both neutralino DM annihilation and pp collisions cannot be the 
mechanism responsible for the HXR excess emission; primary electrons are still a marginally viable solution provided that their 
spectrum has a low-energy cutoff at i5 ^ 30 — 90 MeV. We also find that diffuse radio emission localized at the cluster center 
is expected in all these models and requires quite low values of the average magnetic field (B 0.1 — 0.2 /iC in primary and 
5— i ' secondary-pp models; B ~ 0.055 — 0.39 nG in secondary-DM models) to agree with the observations. Finally, the WR model 
(with B ~ 0.4 — 2.0 /iC) offers, so far, the best description of the cluster in terms of the temperature distribution, heating and 
pressure and multi-frequency spectral energy distribution. Fermi observations of Ophiuchus will set further constraints to this 
model. 
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1. Introduction 

Hard X-Ray (HXR) excess emission in galaxy clusters 
has been observed in the direction of several nearby sys- 
tems (see lNevalainen et al. 2004)) but its origin is still dis- 
puted. It has been proposed that such HXR emission is 
due to inverse Compton scattering (ICS) of relativistic 
electrons with the cosmic microwave background (CMB) 
(Blasi & Colafrancesco 1999, Atoyan & Yolk 2000, Ensslin 
& Biermann 1998, Sarazin 1999, Brunetti 2004, Profumo 
2008, see Petrosian et al. 2008 for a recent review), to 
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bremsstrahlung emission from a supra-thermal electron 
population (Dogiel et al. 2007, see Petrosian et al. 2008 
for a recent review) or to a population of PeV electrons 
that would radiate in hard X-rays through synchrotron 
emission (Timokhin, Aharonian & Neronov 2004; Inoue, 
Aharonian & Sugiyama 2005). None of these models has 
been definitely proven or rejected, so far, due to the lack 
of instrumental sensitivity (spatial and spectral) of the 
available experiments operating in the HXR band. 

In such a context, the Ophiuchus cluster (z=0.028, 
Johnston et al. 1981) has been recently at the center of 
an interesting dispute concerning the combination of new 
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observational evidence for the HXR emission and various 
theoretical considerations on its origin. 
The Ophiuchus cluster seems to have a high plasma 
temperature kT ~ 10 keV (Johnston et al. 1981). 
Measurements of the IC gas temperature vary from 8.5 ± 
0.5 keV (INTEGRAL; Eckert et al. 2008) up to 
keV (Swift/BAT; Ajello et al. 2009). Watanabe et al. 
(2001) also fomid a large (20' x 30'), hot {kT > 13 kcV) 
region, 20 arcmin west of the cluster center, from which 
they concluded that the cluster is not dynamically relaxed, 
and suggested that it experienced a major merging event 
in the recent past {t ^ 1 Gyr). 

Eckert et al. (2008) have recently reported a tentatively 
resolved (^ 5') X-ray source at the cluster center, and in- 
dicated the presence of a non-thermal emission tail with 
a flux (10.1 ± 2.5) • 10-^2 erg cm-^ s"^ in the 20-60 keV 
energy band. These authors interpreted the non-thermal 
hard X-ray emission as due to ICS emission from relativis- 
tic electrons scattered off the CMB in the intra-cluster 
(IC) medium. Suzaku observations of the Ophiuchus clus- 
ter by Fujita et al. (2008) have, however, failed to de- 
tect the non-thermal component detected by Eckert et al. 
(2008), although their quoted upper limit of 2.8- 10""'^'^ erg 
cm~^ s~^ in the 20-60 keV energy band is still compatible 
with the INTEGRAL detection. Ajello et al. (2009) have 
found, using Swift/BAT spectra (with a IC gas tempera- 
ture of kT = 9.5 keV), an upper limit on the Ophiuchus 
non-thermal X-ray emission in the 20-60 keV band, of 
7.2 • 10-12 (.j-g cm-^ (90% c.L). We notice that the 
INTEGRAL detection and the Swift-BAT upper limit are 
consistent, at the same 90% confidence level, in the flux 
range (6.1 — 7.2) • lO^-'^^ erg cm^^ a^^. This is the flux 
range in which the HXR excess detected from Ophiuchus 
is consistent with both Swift-BAT and INTEGRAL ob- 
servations. In our study of the origin of such HXR excess, 
we refer to the value -Pjo-eokeV = 7.2 • lO^^^ erg cm^^ 
s-i as the maximum value of its flux and we discuss how 
our results change by considering also the minimum flux 
of the HXR, -Pjo-eokeV = 6.1 • 10" erg cm-^ s"\ that is 
~ 15% less than the previous maximum value. 

The possible presence of an ICS tail of HXR emis- 
sion in Ophiuchus was previously related to the identifica- 
tion of the steep-spectrum radio source MSH 17-203 (also 
dubbed Cul 1709-231) as a radio mini halo (Johnston et al. 
1981), that would thus imply the presence of relativistic 
electrons, and hence the relative ICS emission emerging 
from the thermal bremsstrahlung emission in the X-ray 
band at ^ > 20 keV. 

However, a recent high resolution study (Perez- Torres et 

al. 2008) made with 240 and 607 MHz GMRT radio ob- 
servations of the Ophiuchus cluster of galaxies, along with 
archival 74 and 1400 MHz VLA data, indicates that there 
is no significant diffuse radio emission in the core of the 
Ophiuchus, and that the previous measurements of radio 
flux from the MSH 17-203 source (Slee & Higgins 1975, 
Slee 1977) do not refer to the radio halo of the cluster: 
these authors present new upper limits to the integrated, 
diffuse non-thermal radio emission from the cluster core. 



More recently, Govoni et al. (2009) pointed out that there 
is indeed a diffuse radio mini-halo located at the center 
of the Ophiuchus cluster, with angular size of ~ 9 x 12 
arcmin^ ; the flux of the mini halo is ~ 8 times lower than 
the old measurement of Johnston et al. (1981), which had 
a resolution of 80 arcmin (see details in Perez- Torres et al. 
2008). 

There is no other information on non-thermal emis- 
sion from Ophiuchus; gamma-ray emission from this clus- 
ter has been not detected and therefore the only infor- 
mation we have directly on the high-E particle popula- 
tion of Ophiuchus is an upper limit obtained by EGRET 
F{> 100 MeV) = 5 • 10"^ cm'^ s'^ (Reimer et al. 2003). 

In such an obervational scenario (that is similar to 
other clusters where an HXR emission detection has been 
claimed) the HXR emission excess from Ophiuchus has 
been recently interpreted as ICS emission from either a 
population of primary cosmic ray electrons (Eckert et 
al. 2008) or secondary electrons produced in neutralino 
DM annihilation (Profumo 2008). In particular, Profumo 
(2008) proposed that a combination of three different 
neutralino DM models [M^ = 81{W+W-),40{bb) and 
10(t"'"t-) GeV] is consistent with all non-thermal emis- 
sion data for Ophiuchus, from radio to HXR and gamma- 
rays. The available data on diffuse radio emission in the 
core of Ophiuchus and the overall analysis of its multi- 
frequency SED further led Perez- Torres et al. (2008) to 
conclude that i) a synchrotron-l-ICS model from primary 
cosmic ray electrons is in marginal agreement with the 
the available data, with a range of magnetic field values 
B ~ 0.02 — 0.3 /iG; ii) that a pure neutralino annihilation 
scenario cannot reproduce both radio and HXR emission, 
unless extremely low magnetic field values (10"^ to 10"'^ 
/iG) are assumed; iii) a scenario in which synchrotron and 
ICS arise from PeV electron-positron pairs (via interac- 
tions with the CMB), can also be ruled out, as it pre- 
dicts a non-thermal soft X-ray emission that largely ex- 
ceeds the thermal bremsstrahlung emission measured by 
INTEGRAL. 

In this paper we take a more radical approach to the 

problem of the HXR emission of Ophiuchus and we con- 
sider not only the SED properties of synchrotron plus ICS 
scenarios (from both primary and secondary electrons) 
but also the physical consequences of the ICS origin of the 
HXR emission in all models so far viable: primary electron 
model (Sect 2.1), secondary electron models from pp colli- 
sions (Sect. 2. 2) and from DM annihilation (Sect. 2. 3) and 
finally a Warming Ray model (Sect. 3). We will discuss our 
conclusions in Sect. 4. 

Throughout the paper, we use a flat, vacuum- 
dominated cosmological model with flm = 0.3, Qa = 0.7 
and h = 0.7. 
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2. Modeling the electron populations in 
Ophiuchus 

The spatial distribution of the intra cluster (IC) thermal 
plasma in the Ophiuchus cluster can be represented by: 



nthir) = nth,o 



1+ - 



-qth 



(1) 



with Tc = lOSh^Q kpc and qth = 0.96 (Watanabe et al. 
2001), and nth,o = 1-77 x IQ-^/i^/^ cm'^ (Johnston et al. 
1981). The available data for this cluster indicate, out to a 
reasonable extent, an isothermal temperature profile with 
a temperature kT ~ 9.9 keV (Watanabe et al. 2001, Ajello 
et al. 2009); the virial radius of the cluster is i? '--^ l.lh^Q 
Mpc (Mohr et al. 1999). 

Non-thermal electrons that can be able to produce the 
cluster HXR emission by ICS on CMB photons must have 
energies Ef, « 0.35 GeV(i?/keV)^/^ in the range « 1.6 — 
2.7 GeV, in the case of the HXR emission observed in the 
20-60 keV range (see, e.g., Colafrancesco et al. 2005). 

In this section we discuss the predictions of various 
models for the origin of the high-E electrons: i) pri- 
mary electron model (PEM) (see, e.g., Colafrancesco, 
Marchegiani & Perola 2005 and references therein) ; ii) sec- 
ondary electron model produced by proton-proton (pp) 
collisions in the cluster atmosphere (SEM-pp) (see, e.g., 
Colafrancesco & Blasi 1998, Blasi & Colafrancesco 1999, 
Marchegiani et al. 2007); iii) secondary electron models 
produced by neutralino Dark Matter annihilation (SEM- 
DM) (see e.g. Colafrancesco et al. 2006, Profumo 2008). 
The electron spectra expected in the previous models are 
normalized by assuming that the produced ICS HXR emis- 
sion equals the Swift-B AT/INTEGRAL data. 

In addition to the previous models, we also consider in 
the next Sect. 3 a self-consistent warming-ray (WR) elec- 
tron model (see e.g Colafrancesco & Marchegiani 2008; 
Colafrancesco et al. 2004) in which the cluster atmosphere 
is heated - in a quasi stationary equilibrium condition be- 
tween heating and cooling - by the interactions of non- 
thermal cosmic-ray protons with the IC gas. Note that 
this WR model reproduces the X-ray properties of the 
thermal IC gas (namely its temperature and density pro- 
files) and, therefore, we use this constraint to predict the 
cluster ICS HXR emission. 



2.1. Primary Electron Model (PEM) 

The high-E electron spectrum is usually best constrained 
by using the radio halo synchrotron spectrum that pro- 
vides direct information on the electron spectral shape (see 
e.g. Colafrancesco et al. 2005 for a discussion). However, 
the upper limits on the Ophiuchus radio halo (see Perez- 
Torres et al. 2008) and the measurement at 1.4 GHz made 
by Govoni et al. (2009) are not sufficient to determine 
precisely the electron spectrum. Therefore, we choose to 
adopt here a simple power-law model 



with p = 3.5, which corresponds to a radio spectral index 
of an = 1.25, typical of radio halos in galaxy clusters and 
consistent with the available radio data on Ophiuchus; 
we also consider the case with p = 4.4, which provides 
an = 1.7, which is, at la level, the maximum value of the 
spectral index allowed by radio data (see, e.g., FigH]). The 
radial function 17(7') should be constrained, in principle, by 
measuring the radial shape of the ICS emission of the rel- 
ativistic electrons (the radio measures are not sufficient to 
constrain this shape; see discussion in Colafrancesco et al. 
2005); since the available HXR measures of the Ophiuchus 
cluster have no sufficient spatial resolution, we assume in 
the following that the radial distribution function of the 
non-thermal electrons has the same shape of the thermal 
gas radial distribution (see Eq[T]). Under such assumption, 
we consistently assume that the relativistic electrons ex- 
tend out to the virial radius of the cluster, like the IC 
plasma. 

The value of the normalization of the electron spectrum 
Neo can be derived by reproducing the value of the 
HXR flux set by the Swift-BAT and INTEGRAL exper- 
iments: for the previous spectral index, we obtain the 
value A^e,o = 1-1 x 10"^" GeV-^ cm'^ for p = 3.5, and 
Ne,o = 2.1 X 10-^° GeV-1 cm^^ for p ^ 4.4. 

All the previous information can be derived, strictly 
speaking, only for the range of the electron energies that 
produce the HXR emission via ICS. To obtain information 
on other energy ranges of the electron spectrum one must 
consider other constraints. 

The only other constraint that can be set on the elec- 
tron spectrum comes from the requirement that the heat- 
ing rate of the IC gas produced by non-thermal electron 
Coulomb collisions does not exceed the bremsstrahlung 
cooling rate of the IC gas. The heating rate produced by 
an electron with Lorentz factor 7 and velocity v — (3c is 
given by 



dE 



1 



In 



2 me 



(3) 



where is the (suitably averaged) squared charge of the 



plasma's nuclei, K — Att nth ft Tn^ 



2.82 fm, and Ir, ^ he 
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2004, 



the plasma frequency (see Colafrancesco et al. 
Colafrancesco & Marchegiani 2008). 

As a consequence, the heating rate induced by the elec- 
trons with the spectrum assumed in Eq.([2]) is given by: 



rfe 



dE 



(4) 



while the cooling rate is given by 

de 
dt 



a[ntbir)f VkT{r, t); 



33 h ml 



Ne{E,r)^Ne,o{E/GeY)-P-g(r) 



(2) 



4.8 X lO^^^'z- 



1 erg cm'^ 



(5) 
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where z is an average charge of the IC plasma (we have 
approximated here the Gaunt factor G by unity). 

The two expressions in Eqs. ^ and (O equal (to the 
value dtjdt ^ 4.7 x 10^^^ erg cm^"^ s^^) for the energy 
^ 33 MeV for p = 3.5, and for £;,„i„ - 91 MeV 
for p = 4.4; therefore, we assume these energy values as 
the minimum energy of the primary electron spectra. As 
for the maximum electron energy we can safely choose 
Emax — >■ oo since its specific value is irrelevant for the as- 
sumed spectral index. 

Under such assumptions, it is possible to calculate the 
overall radiation emission via the various emission mech- 
anisms in different frequency ranges. 

The synchrotron emission spectrum produced at radio 
frequencies by primary electrons is shown in Fig [T] for dif- 
ferent values of the magnetic field (the B-field has been 
assumed to be constant in the emission region, and this 
corresponds, approximately, to consider a volume aver- 
aged value of the magnetic field). From the radio emission, 
we can derive a value of the average magnetic field of ~ 
0.1 fiG for p = 3.5 and ~ 0.2 fj,G for p = 4.4, in agreement 
with that derived by Ajello et al. (2009) and Perez-Torres 
et al. (2008). 

In the gamma-ray frequency range, these primary 
electrons emit via non-thermal bremsstrahlung and ICS 
against CMB radiation field if their energy spectrum, 
as we assume in this case, extends up to high energies 
(at least up to 100 — 1000 GeV in order to pro- 

duce ICS emission in the energy range 0.1-10 GeV). Fig. 
[2] shows that the EGRET upper limit on Ophiuchus, 
F{> 100 MeV) < 5 X 10"* cm^^ g-^ (Reimer et al. 2003), 
is not exceeded in the p ~ 3.5 case, while it is marginally 
exceeded in the p = 4.4 case. In the first case, we can con- 
clude that the HXR observation of Ophiuchus cluster sets 
a constraint on the ICS emission from relativistic electrons 
that is stronger than the analogous limit set by EGRET; 
in the second case, the EGRET limit is stronger than the 
HXR limit. The signals we derive here for the gamma-ray 
emission of Ophiuchus in the PEM, and in particular the 
one derived from non-thermal bremsstrahlung emission, 
are sensibly larger than the Fermi sensitivity at E ^ 300 
MeV; therefore, such an experiment could be able either 
to detect the bremsstrahlung gamma-ray emission from 
Ophiuchus or set even stronger limits on the non-thermal 
electron density. 

It is important to stress here that the HXR data sets 
also indirectly a lower cut-off of E„iin ~ 33 and 90 MeV 
for the two models considered on the electron spectrum in 
order to have an heating rate not larger than the cooling 
rate. 
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Fig. 1. The diffuse radio emission spectrum as produced from 
primary electrons with p — 3.5 (upper panel) and p = 4.4 
(lower panel) is shown for different values of the uniform mag- 
netic field (in units of /iC), as labelled. Data are from Perez- 
Torres et al. (2008) (upper limits) and Govoni et al. (2009) 
(point at 1.4 GHz). 



(see Marchegiani et al. 2007 for details), we assume that 
the non-thermal protons have the following spectrum 



NpiE,r)^Np,o{E/GeYr'g{r) 



(6) 



2.2. Secondary Electron Model from pp collisions 
(SEM-pp) 

To calculate the overall radiation emission from secondary 
electrons produced by collisions of cosmic ray (CR) non- 
thermal protons and the thermal protons of the IC gas 



and we further assume, also in this case, that their spatial 
distribution is the same of the thermal IC gas out to the 
virial radius. 

Similarly to the PEM model, we assume proton spectral 
indices s = 2.5 and s — 3.4, which provide again radio 
spectral indices of = 1.25 and 1.7 respectively. 
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Fig. 2. The gamma- ray spectrum produced by primary elec- 
trons with p = 3.5 (upper panel) and p = 4.4 (lower panel) 
via ICS (dashed) and bremsstrahlung (dot-dashed) are com- 
pared to the sensitivity curves of EGRET and Fermi {5a, 1 
year observation). 

The HXR data explained in terms of the ICS emis- 
sion of the secondary electrons, provides values iVp.o = 
l.SxlO-^GeV-i cm-3 (for s = 2.5) and iVp,o = 3.4x10-5 
GeV^^ cm-'^ (for s — 3.4) for the normalization of the 
spectrum in Eq.Q. 

Such values of the proton density normalization imply se- 
rious problems for Ophiuchus: 

i) first, the proton pressure at the center of the cluster is 
more than 5 times larger than the thermal gas pressure 
for s — 2.5 and 367 times larger for s = 3.4, a fact that 
sets serious problems to the cluster stability; 

ii) secondly, the heating rate (Coulomb losses and 



hadronic collisions) induced by non-thermal protons in 
the cluster center is de/dt ^ 1.1 x 10^^^ erg cm^'' s^^ 
for s ~ 2.5, and de/dt ~ 1.4 x 10^^"' erg cm~^ s^^ for 
s — 3.4; these values are about 23 and 298 times larger 
than the cooling rate. This fact would imply a quite fast 
heating of the cluster, d[kT)/dt ~ 41 keV Gyr~^ and 
d{kT)/dt - 531 keV Gyr-\ that will bring in a short 
time the whole IC gas to a temperature sensitively differ- 
ent (larger) from the observed one; 

iii) finally, the gamma-ray emission produced by both sec- 
ondary electrons and by neutral pion decay (see Fig[3]) 
exceeds the EGRET upper limit on Ophiuchus by a fac- 
tor ~ 18 and 170, for s = 2.5 and 3.4 respectively. 

Thus, we must conclude that the HXR emission of 
Ophiuchus as set by Swift and INTEGRAL, cannot be 
produced by secondary SEM-pp electrons. We notice that 
such a conclusion is analogous to that found in the case 
of other clusters we already studied like Coma, A2199, 
A2163 and Perseus (Colafrancesco & Marchegiani 2008). 

2.3. Secondary Electron Model from DM annihilation 
(SEM-DM) 

We consider here three neutralino DM models (similarly 
to the analysis of Profumo 2008) with neutralino masses 
= 81(iy+iy-),40(&6) and 10(r+r-) GeV. 
For each neutralino model we consider a radial DM 
density profile as given by 

goMir) = exp[~(2/a)((r/r,)" - 1)] (7) 

(Navarro et al. 2004), with a = 0.17 and Tc equal to the 
core radius of the thermal gas density distribution. We as- 
sume that this DM radial profile extends out to the virial 
radius. The spectrum of the DM source function for the 
secondary electrons has, consequently, a radial distribu- 
tion cx 5DAfW- 

To derive the equilibrium spectrum of these secondary 
electrons in Ophiuchus we consider the role of the domi- 
nant energy loss mechanisms. These are ICS losses against 
CMB photons and synchrotron losses for electrons with 
energy larger than a few hundreds MeV (notice that syn- 
chrotron losses for magnetic fields less than 3 /iG, are neg- 
ligible with respect to the ICS losses), while at low energies 
( ^ 150 MeV) the dominant energy loss mechanisms are 
Coulombian interactions with the IC gas particles. 
For this reason the final spatial distribution of secondary 
electrons is proportional to high energies (> 

150 MeV) and proportional to sIjm ('') /'^t'» low en- 

ergies (< 150 MeV). 

The DM-produced secondary electron density is fixed, 
also in this case, by requiring that their ICS emission fits 
the observed HXR emission; such a constraint corresponds 
to set the value of the neutralino annihilation cross section, 
(crF), because both the neutralino mass and its composi- 
tion have been fixed by the chosen model. 

Also this SEM-DM model has serious implications for 
the Ophiuchus cluster. 
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Fig. 3. The gamma-ray spectrum of Ophiuchus with s — 2.5 
(upper panel) and s = 3.4 (lower panel) as produced by ICS 
(dashes) and bremsstrahlung (dot dashes) of secondary elec- 
trons and by neutral pion decay (long dashes). We compare the 
predictions of the SEM-pp model with the sensitivity curves of 
EGRET and Fermi {5a, 1 year observation). 



i) The heating rate at the cluster center as produced by the 
secondary SEM-DM electrons is very high; Fig. [4] shows 
the secondary electrons heating rate at different radii com- 
pared to the cooling rate of the thermal IC gas. In fact, 
the heating rate largely exceeds the cooling rate in the 
cluster core at r < 30 kpc. This result would imply a 
fast over-heating of the Ophiuchus core, even though the 
volume integral of the heating rate is always lower than 
the volume integral of the cooling rate (this last quan- 
tity is 1.5 X 10^^ erg/s, while the integrate heating rate is 



4.6 X 10^'^ crg/s, 7.4 x lO*^ erg/s and 3.6 x lO''^ erg/s for 
= 81, 40 and 10 GeV, respectively). 
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Fig. 4. The heating rate induced by secondary electrons pro- 
duced by DM annihilation is shown at different radii for dif- 
ferent neutralino models: — 81 GeV (dashed), 40 GeV 
(dot-dashed) and 10 GeV (dot-dot dashed). The intracluster 
gas cooling rate (solid) is also shown for comparison. 

ii) We show in Fig. [5] the gamma-ray emission spec- 
tra as produced by the DM composite model worked out 
here via the three main mechanisms of gamma-ray emis- 
sion: ICS and bremsstrahlung from secondary electrons 
and neutral pion decay. All the three DM models con- 
sidered in this composite DM model for Ophiuchus pro- 
duce a gamma-ray flux that exceeds the EGRET limit, 
F(> 100 MeV) = 5.0 x 10"^ cm^^ g-^ (the low mass 
neutralino model with = 10 GeV is marginally con- 
sistent with the EGRET limit). The gamma-ray flux of 
Ophiuchus at > 100 MeV produced under the assump- 
tion that the same DM model reproduce the HXR data 
are 7.4 x 10~^, 1.3 x 10"'^ and 4.3 x 10'^ cm"^ s^^ for neu- 
tralino masses of 81, 40 and 10 GeV. A direct prediction 
of this DM model is that the gamma-ray flux produced 
by the three neutralino models considered here should be 
easily detectable by the Fermi experiment whose results 
will be able, therefore, to validate or rule out this model 
for the origin of the HXR emission of Ophiuchus. 

iii) Fig. [S] shows the diffuse synchrotron radio spec- 
tra produced by the same secondary SEM-DM electrons 
under the assumption of a reference value of the average 
intracluster magnetic field in Ophiuchus of 0.1 /iG. This 
figure shows that, for this value of magnetic field, a model 
whith neutralino mass between 40 and 80 GeV can repro- 
duce the radio data. 

We also searched for the value of the magnetic field that, 
for each of the models we considered, reproduces the 
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Fig. 5. The overall gamma-ray spectrum produced by the 
composition of the contributions of the secondary SEM-DM 
electrons (ICS and bremsstrahlung) and of the neutral pion 
decay for the three neutralino models here considered. The sen- 
sitivity curves of EGRET and Fermi (5a, 1 year observation) 
are shown for comparison. 



Fig. 6. The radio halo spectrum produced from secondary 
SEM-DM electrons via synchrotron emission in a constant 
magnetic field of 0.1 piG. Data are from Perez- Torres et al. 
(2008) (upper limits) and Govoni et al. (2009) (point at 1.4 
GHz). 



Ophiuchus radio halo flux at 1.4 GHz. Fig[7] shows the 
radio spectrum produced by SEM-DM electrons for best- 
fit magnetic field values of 0.055, 0.18 and 0.39 //G, for 
= 81, 40 and 10 GeV, respectively. We can conclude 
that the 81 GeV (W^+W^~) model is consistent with radio 
data, while the 10 GeV (t+t^) model is not consistent. 
The 40 GeV (bb) model is a border-line situation: we find, 
in fact, that for a slightly lower magnetic field of 0.17 /iG 
the radio spectrum is marginally consistent with the point 
at 1.4 GHz and the upper limit at 74 MHz. 

3. Warming Ray Model 

In this Section we abandon the strategy of fitting the ICS 
emission produced by high-E electrons to the HXR data of 
Swift and INTEGRAL because we evaluate - in the frame- 
work of the Warming Ray (WR) model (see Colafrancesco 
& Marchegiani 2008, Colafrancesco et al. 2004) ~ the spec- 
tral and spatial characteristics of the WR proton popula- 
tion that produce through their heating action the tem- 
perature structure of the Ophiuchus cluster, namely a con- 
stant temperature radial profile at the observed value of 
kT « 9.9 keV. 

The proton spectrum is written as in Eq. with 
the values s — 2.5 and s = 3.4, and assuming a radial 
distribution given by g(r) cx g"j^{r), where the value of a 
is found by fitting the radial temperature profile of the 
cluster (see Colafrancesco & Marchegiani 2008 for techni- 
cal details). The best fit analysis of the temperature pro- 
file of Ophiuchus provides the value a = 1, in analogy 



to what is found for other isothermal clusters (see dis- 
cussion in Colafrancesco & Marchegiani 2008), while the 
central WR density is Np^Q = 4.9 x 10^* cm^^ GeV^^ 
and A^p,o = 9.4 x 10"^ cm'^ GeV~i for s = 2.5 and 
s — 3.4 respectively, i.e. a factor « 31 and 362 lower 
than that required to reproduce, in this model, the Swift 
BAT/INTEGRAL HXR data. Consequently, the pressure 
ratio of the WR to the thermal gas at the cluster center 
is PcRlPth ~ 0.17 for s = 2.5 and PcrIPik ~ 1-0 for 
s = 3.4: the first value does not give any problem to the 
overall stability of the cluster, while the second one is a 
problematic situation. 

The diffuse radio emission produced by the secondary 
electron in this WR model is shown in Fig. [5] for various 
values of the uniform magnetic field; our results indicate 
that a uniform B-field of the order of ^ 0.4 fiG is required 
to fit the available data for s = 2.5, and B ~ 2.0 /iG for 
s = 3.4. 

Fig. [9] shows the diffuse gamma-ray emission from 
Ophiuchus as expected in the WR model: this emis- 
sion consists of the combination of the neutral pion de- 
cay spectrum and ICS and bremsstrahlung emission from 
secondary electrons. The overall gamma-ray emission of 
Ophiuchus in the WR model, F{> 100 MeV) = 3.3 x 10"^ 
cm~^ s~^ and 2.3 x 10~^ cm"^ s^^ for s = 2.5 and s = 3.4 
respectively, is below the EGRET limit. However, the neu- 
tral pion decay gamma-ray emission predicted in this WR 
model should be detectable by Fermi in 1 yr observation 
(5cr C.I.). 
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Fig. 7. The radio halo spectrum produced from secondary 
SEM-DM electrons via synchrotron emission in a constant 
magnetic field of 0.055, 0.18 and 0.30 //G for Af^ = 81, 40 
and 10 GeV, respectively. Data are from Perez-Torres et al. 

(2008) (upper limits) and Govoni et al. (2009) (point at 1.4 
GHz). 

4. Discussion and conclusions 

We have show in this paper that various and serious prob- 
lems stand out with the ICS interpretation of the HXR 
excess emission of the Ophiuchus cluster and, in general, 
of galaxy clusters for which an HXR emission excess has 
been detected. These problems are: 

i) the actual level of the HXR emission flux: the deriva- 
tion of an HXR emission excess in clusters seems to depend 
strongly on the precise determination of the background 
thermal bremsstrahlung emission. Eckert et al. (2008) de- 
rived from INTEGRAL data an HXR flux in the 20-60 
keV band of Frxr = (10.1 ± 2.5) x lO^^^ erg cm-^ s-\ 
with an IC gas temperature of S.SG^Q gg keV. Ajello ct al. 

(2009) , using Swift BAT data, derived an upper limit (90% 
c.l.) Fhxr < 7.2 X 10'i2 erg cm'^ s'S with a different IC 
gas temperature of 9.5l}'^ keV; it must be noticed that 
the same authors (Ajello et al. 2009) derived the upper 
limit for the HXR flux of Frxr < 4.5 x IQ-^^ erg cm^^ 
s""'^ by using an higher value of the temperature 9.931q'24 
keV, as obtained by using a combination of Chandra and 
Swift-BAT data. 

In conclusion, it is clear that a crucial input quantity to 
determine the value of the HXR excess flux is the detailed 
modeling of the thermal emission of the IC gas, because 
different values assumed for the IC gas tempcratiirc lead 
to different conclusions on the amount of the HXR excess 
flux (see, e.g. the long standing discussion on the evidence 
and counter evidence of the HXR excess in Coma, Fusco- 
Femiano et al. 1999, 2004, 2007, 2008; Rossetti & Molendi 
2004, 2007; see also Petrosian et al. 2008 for a review). For 




0.001 I , , , , 

10 100 1000 10000 
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Fig. 8. The radio halo spectrum produced by secondary elec- 
trons with s = 2.5 (upper panel) and s = 3.4 (lower panel) in 
the WR model is shown for different values of a constant IC 
magnetic field (as labeled) in units of /uG. Data are from Perez- 
Torres et al. (2008) (upper limits) and Govoni et al. (2009) 
(point at 1.4 GHz). 



this reason, it would be extremely important to estimate 

the temperature of the IC gas through measurements that 
are independent from those obtained in the X-ray band. 
We notice, in this context, that a reliable method to mea- 
sure IC gas temperatures can be found by using Sunyaev- 
Zel'dovich observations over a wide spectral band (from 
^ 100 to ^ 400 GHz) reaching high frequencies where 
the sensitivity of the SZE to the cluster temperature is 
maximum (we have discussed in details this issue in a 
dedicated paper, Colafrancesco & Marchegiani 2009; see 
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Fig. 9. The gamma-ray spectrum produced by secondary elec- 
trons with s — 2.5 (upper panel) and s — 3.4 (lower panel) 
via ICS (dashes) and bremsstrahlung (dot-dashes) and neutral 
pion decay (solid) in the WR model that is set up by the re- 
quirement to reproduce the observed temperature structures 
of the IC gas in Ophiuchus. The sensitivity curves of EGRET 
and Fermi (5(t, 1 year observation) are shown for comparison. 

also Colafrancesco et al. 2003 and Colafrancesco 2007 for 
a review); 

ii) the ICS HXR scenario: the hypothesis that high- 
E electrons are responsible for an ICS HXR emission at 
the level found by the combination of Swift-BAT and 
INTEGRAL observations leads to some important con- 
sequences. 

First, in order to reconcile the HXR excess value with 
the relative diffuse synchrotron radio emission (from the 
same electron population) at the level observed in the 



same cluster, the value of the average magnetic field must 
be quite low and of the order of ~ 0.1 for p = 3.5 
and ^ 0.2 fiG for p = 4.4 (see Fig. [T]). The result 
found for Ophiuchus is analogous to what is derived for 
other clusters (see also our previous results discussed in 
Colafrancesco, Marchegiani & Perola 2005, Marchegiani, 
Perola & Colafrancesco 2007): Specifically, we found that 
the data are consistent with an IC magnetic field of order 
of ~ 0.7 and 1.2 /iC at the cluster center with a decreasing 
radial profile similar to that of the IC gas, for p = 3.5 and 
4.4 respectively. 

Secondly, there is a strong relation between the ICS HXR 
emission level and the relative gamma-ray emission and 
the consequences on the physics of the cluster, i.e. the 
heating of the IC gas and the ratio between non-thermal 
and thermal pressures: 

- if the electrons that produce the HXR emission are pri- 
maries, their gamma-ray emission (dominated at i? < 1 
GeV by non-thermal bremsstrahlung) is slightly lower 
than the EGRET upper limit in the p = 3.5 model, 
and slightly higher than this limit in the p — 4.4 model, 
but certainly detectable by Fermi (see Figl2]). If Fermi 
will not detect such gamma-ray emission, one should con- 
clude that the ICS HXR emission is much lower than the 
Swift- INTEGRAL HXR detection and that the relativis- 
tic electron content of Ophiuchus is consequently much 
lower. The HXR data set, in addition, a lower cut-off of 
Emin ~ 33 and 93 MeV (for p =3.5 and 4.4) on the elec- 
tron spectrum in order to have an heating rate not larger 
than the cooling rate; 

- if the electrons responsible for the ICS HXR emission are 
secondary particles produced in the decay of charged pi- 
ous generated by cosmic-ray proton collisions with the IC 
gas protons (SEM-pp), then an ICS flux set at the HXR 
observations has unacceptable consequences. Specifically 
we find that: in the s = 2.5 case, the pressure exerted 
by relativistic protons at the cluster center is ~ 5 times 
larger than the thermal gas one; the heating rate induced 
by the same relativistic protons at the cluster center is 
~ 23 times larger than the IC gas cooling rate; and the 
gamma-ray emission produced by neutral pion decay ex- 
ceeds the EGRET limit by a factor ~ 18; in the s = 3.4 
case, these quantities rise respectively to ~ 367, 298 and 
170 (see FiglS]). For all these reasons, we conclude that if 
electrons produce an ICS HXR emission in the observed 
range, they cannot be secondary (in the SEM-pp). This 
conclusion is analogous to what has been found also in 
other clusters like Coma, A2199, A2163 and Perseus (see 
Colafrancesco & Marchegiani 2008); 

- if electrons are produced by neutralino DM annihila- 
tion, we have found that: the heating rate they induce 
at the cluster center is quite high (see Fig.H]); the relative 
gamma-ray emission exceeds the EGRET limit for the two 
high-Af^ models here considered {My, — 40 and 81 GeV) 
(see Fig. [5]) , with a marginal consistency for the low-mass 
model with M^, = 10 GeV; the radio flux produced by 
electrons is consistent with the available data for M^^ > 40 
GeV, and for B < 0.18 /LtC. Therefore, the information in- 
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ferred by gamma and radio data are not compatible, and 
we conclude that it is not possible to conceive that the 
ICS HXR emission of secondary SEM-DM electrons has a 
flux close to the available observation (i.e., the maximum 
allowed flux set by Swift and INTEGRAL, see Sect.l), and 
thus their annihilation cross section must be much lower 
than the values used by Profumo (2008). 
Even normalizing these models to the lower allowed flux 
value of the HXR excess of Ophiuchus (see our discus- 
sion in Sect.l), all the previous results vary (decrease) by 
~ 15%, leaving unchanged our basic conclusions. 

iii) Relaxing the assumption to recover the observed 
HXR excess and assuming that non-thermal protons act 
as warming rays (see Colafrancesco & Marchegiani 2008) 
it is possible to paint a much more acceptable scenario 
in which the unacceptable pressure ratios derived in SEM 
models do not hold since the ratio Pnon-th/ Pth ~ 0.17 and 
1.0 for, respectively, s — 2.5 and 3.4 and it is constant 
throughout all the cluster (this is because non-thermal 
protons must have the same spatial distribution of the 
thermal IC gas to recover the spatial temperature distri- 
bution of the cluster). In addition, the WR model has 
other positive aspects for the cluster structure: i) it does 
not induce excess heating effects, since a quasi-stationary 
balance between heating and cooling is the working as- 
sumption of the WR model; ii) we found that the diffuse 
radio emission produced in this case requires, for s = 2.5 
and s — 3.4 respectively, a value of the average magnetic 
field of ~ 0.4 and 2 /xG (see Fig. |8]) and a central value of 
'-^ 1.1 and 6 /iG with a radial profile similar to that of the 
IC gas, consistently with the general findings for clusters 
through Faraday Rotation measurements (see, e.g., Carilli 
k Taylor 2002, Govoni & Feretti 2004); ui) the gamma- 
ray emission produced in this model is quite lower than 
the EGRET limit but definitely detectable by Fermi (see 
Fig. [9]). The Fermi detection of such gamma-ray emission 
from Ophiuchus will have a crucial impact for proving or 
disproving this model. 

In such a WR model, the HXR ICS fiux of Ophiuchus is 
much lower (by a factor ^ 30 and 362 for s = 2.5 and 
s = 3.4 respectively) than the limit set by the present 
observations (by INTEGRAL and Swift-BAT) and could 
only be detectable by using long exposure observations 
with the next generation HXR instruments like NeXT (see 
e.g. Takahashi et al. 2004). 

To conclude, models of high-E electrons in clusters that 
can be adjusted to reproduce their ICS HXR emission 
at the level indicated by the available observations fail 
to work because they would imply unacceptable levels of 
heating and gamma-ray emission. On the contrary, models 
of high-energy particles that are able to reproduce the IC 
gas temperature structure (i.e. that WR model) predict a 
level of non-thermal HXR ICS emission that is far below 
the current limits obtained with INTEGRAL and Swift- 
BAT, and provide an overall Spectral Energy Distribution 
that is consistent with all the available data ~ from radio 
to gamma-rays - on Ophiuchus as well as on other clusters. 
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